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ABSTRACT 

Star clusters arc formed in molecular clouds whi ch are believed to be the birth places 
of most stars. From recent observational data, iLada fc Ladal (|2003( ) estimated that 
only 4% to 7% of the clusters embedded inside molecular clouds have survived. An 
important mechanism for the disruption of embedded (bound)-clusters is the disper- 
sion of the parent cloud by UV radiation, stellar winds and/or supernova explosions. 
In this work we study the effect of this mechanism by N-body simulations. We find 
that most embedded-clusters survive for more than 30 Myr even when different initial 
conditions of the cluster may introduce some minor variations, but the general result 
is rather robust. 

Key words: stellar dynamics - methods: N-body simulations - galaxies: star clusters. 



1 INTRODUCTION 

Stellar clusters are among the most interesting objects in 
astronomy. We are interested in how they form, how they 
evolve and how they die. Observations gave the evidence 
that most stars are not born independently but in stellar 
clusters or stell ar associations which are formed in molecu- 
lar clouds (e.g., iLada fc Lad"al (|2003l )). Due to the limits of 
observational techniques, we do not know in detail the rela- 
tion between the clusters and the parent molecular clouds. 
It is believed that stars are born in clusters and become field 
stars after the clusters disassociate. 

Recently, near infrared observational data (2MASS, 
Two Micron All Sky Survey project at IPAC/Caltech) have 
shown that the number of embedded clusters is much higher 
than the number of optical clusters for which the parent 
clouds have already dissipated and that the survival proba- 
bility for embe d ded c lusters in Milky Way is about 4% to 7% 
(|Lada fc Ladal (|2003l )). Similar evidences for infant mortal- 
ity are found in Ant ennae galaxies (iFall et al.l (|2005l ')'). Small 
Magellanic Cloud (IChandar et alj l|2006l )') and NGC1313 
IPellerin et al.l (|200m 



These imply that clusters are likely to be disrupted be- 
fore the clouds are dissipated completely. Since the time the 
clusters were born, they are under constant threats from 
their surrounding environment. 

Galactic tidal forces , close encou n ters w ith giant molec- 
ular clouds (see, e.g., ICieles et al.l (|2006l ')). shock heat- 



ing and mass l o ss by massive member stars (see e.g., 
iBoilv fc Kroupal (|2003al lbh) are possible disruption mech- 
anisms. Nonetheless, most of these mechanisms have a 
destruction timescale longer than the upper limit of the 
lifetime of molecul ar clouds which is about few to few 
tens Myr (see e . g.. iBUtz fc Shul Il98^ lElmegreeiil (|2000D . 
iHartmann et all (|200ll ). iBonnell et all (|2006l ')'). The results 
of CO observations in the Galaxy suggest that the lifetime 
of molecular clouds i s of the order of 10 Myr (see, e.g., 
iLeisawitz et al.l (|l989l )). The estimation of the timescale for 
photoevaporation and statistics on the expected numbers of 
stars per cloud show that giant molecular clouds of mass 
10^ Mp) are expected to su rvive for about 30 Myr (see, e.g., 
IWiUiams fc McKed (|l997h ). 

Generally speaking, mechanisms with destruction 
timescales less than the lifetime of the clouds should be re- 
sponsible for the low survival probability mentioned above. 
In this work, we focus on the role played by the dispersion 
of the parent cloud on the early evolution of the embedded 
clusters. 

In the beginning, the clusters are bound to their parent 
molecular cloud. As the cloud dissipates, the binding energy 
from the cloud decreases and the stellar systems become out 
of equilibrium. Once out of equilibrium, they may expand 
or dissociate completely. 

The survivability of clusters under gas dispersion 
has been examined in the past and extensive N-body 
sim ulations we r e per fo rmed in t he pa st few years (s ee 
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e.g Lada ot "a ll (119841). ICoodwinl Il997l), iBoily fc Kroupal 
2003 b). Baui ngardt fc Kroupal (|2007| ). ISastian fc GoodwinI 



200d ). iGoodwin fc BastianI (|2006h V In a large set of sim- 
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ulations, iBaumgardt fc Kroupal (|2007l ) studied the disper- 
sal of the residual gas by decreasing the mass with differ- 
ent star formation efficiency (SFE), and in different tidal 
fields. They concluded that the clusters had to form with 
SFE ^ 30% in order to survive gas expulsion, and the ex- 
ternal tidal fields have significant influences only if the ra- 
tio of half mass rad i us to tidal radius is larger than 0.05 . 
iGoodwin fc BastianI (120061 ) and lBastian fc GoodwinI (120061 ) 
addressed similar problem and found that the embedded 
clusters would be destroyed within a few tens of Myr if the 
"effective star formation efficiency" < 30% (eSFE= 1/2Q, 
and Q=0.5 for virial equilibrium). 

The paper is organized as follows. In Section 2, we de- 
scribe the model and parameters for our simulations. In Sec- 
tion 3, we present and discuss the results and statistics of 
the simulations. A summary and some remarks are provided 
in Section 4. 



2 MODEL AND SIMULATION 

We intend to learn the behaviour of a star cluster in a dis- 
persing molecular cloud. For simplicity, we do not consider 
the feedback of the cluster onto the cloud, and the cloud is 
simply represented by its gravitational force. In other words, 
we study the behaviour of a star cluster in a time varying 
gravitational field. We adopted the N- body simulation code 
NB0DY2 developed bv lAarsethI (|200lh for our calculations. 
As a first attempt, we assume a spherically symmetric ex- 
ternal potential (representing the molecular cloud) , and ini- 
tially the centre of mass of the cluster coincides with the 
centre of the potential. The initial spatial and velocity dis- 
tributions of the cluster depend on the initial profile of the 
external potential (i.e., the initial mass and compactness of 
the cloud). Subsequent evolution of the cluster depends, of 
course, on its initial distribution and the rate of dispersion 
of the external potential. 



2.1 Model for the cloud 

The cloud is represented by a Plummer potential, 
-Ghh 



$p — 



(1) 



where Mh is the total mass of the cloud, a is the length 
scale of the potential and G is the gravitational constant. To 
model the dispersion of the cloud we consider the potential 
to evolve in time according to 



■ doe 



(2) 



where ao and a are the initial length scale and the dispersion 
rate of the cloud, respectively. The total mass of the cloud 
remains constant as the length scale increases with time. 
Stellar masses would change by stellar evolution. However, 
since we run our simulations for 30 Myr only (which is about 
the maximum life time of molecular clouds) and massive 
stars are rare, hence we do not consider stellar evolution in 
this work. 

There are three parameters for the cloud; (i) dispersion 
rate a, (ii) total mass Mh, and (iii) initial length scale oq. 



Table 1. Dispersion rates of the cloud 

index A B C D 



OL 

te [Myr] 



0.1 
3.3 



0.2 
1.5 



0.3 
1.1 



0.4 
0.75 



0.5 
0.625 



Mry 



1 Mry 



3 t^r 



2 



r [pc] 
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r [pc] 



5 Mry : 
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Figure 1. Evolution of the cloud potential in dispersion rate of 
a=0.1(dotted line), 0.3 (dashed line) and 0.5 (dash-dotted line), 
whose e-fold times are 4^=3.3, 1.1 and 0.625 Myr, respectively. 
At Myr, the potentials are the same. After 1 Myr, every one 
becomes shallower in difference degree. At 3 Myr, the two cases 
with te=l.l and 0.625 Myr approach zero. At 5 Myr, even the 
case with te=3.3 Myr becomes rather flat. 



2.1.1 Dispersion rate a 

We take a as 0.1, 0.2, 0.3, 0.4, 0.5 in the unit used by the 
code. These correspond to a e-fold time te as 3.3, 1.5, 1.1, 
0.75, 0.625 Myr in real unit (see Table [l|. 

Fig. [T] shows how the cloud potential evolves. The po- 
tentials of a = 0.3 and 0.5 are almost zero after 5 Myr. 
Even for a = 0.1, the potential is rather flat after 5 Myr. 
We would expect the cloud exerts no effect on the cluster 
after a relatively short time in these dispersion rates. The 
seed of destruction is planted (if at all) only in the early 
dispersion stage of the cloud. 



2.1.2 Total mass A4t 

The more massive the cloud is the more influence on the 
cluster is. We consider the total mass of the cloud Mi, from 
0.5 to 10 cluster mass Mc as listed in Table [21 We note that 
Afc ~ 25OOM0 in all our simulations. Generally, the star 
formation efficiency is defined as 

"^m-TmI- 

The mass range we choose corresponds to 77 = 50 to 9%. 
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Table 2. Total cloud mass in unit of cluster mass 



index 


01 


02 


03 


04 


05 


06 


07 


Mb [MJ 


0.5 


1 


1.5 


2 


2.5 


3 


3.5 




08 


09 


10 


11 


12 


13 


14 


Mt [MJ 


4 


4.5 


5 


5.5 


6 


6.5 


7 




15 


16 


17 


18 


19 


20 




Mf, [MJ 


7.5 


8 


8.5 


9 


9.5 


10 





Table 3. Initial length scale of the cloud 



index 


a 


b 


c 


d 


e 


f 


g 


ao [pc] 


0.125 


0.25 


0.375 
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0.625 


0.75 


0.875 
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i 


j 


k 


1 


m 


n 


ao [pc] 
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1.125 
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1.375 
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1.625 


1.75 
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Figure 2. Contour of /3, the cluster-cloud mass ratio at r4o of 
the clusters. Cloud and cluster have the same mass within r4o 
when (i = 0.5. 



For later discussions, we introduce the cluster-cloud 
mass ratio at 40% Lagrangian radius of the cluster 



Mc,r40 + Mb. 



(4) 



where Mb,r4,o is the mass of the cloud within the initial 40% 
Lagrangian radius of the cluster. Mb^rio depends on ao and 
Mb- Fig. [5] shows the contour of /9 in the parameter space 
{Mb,ao). Note that Mc, r4o is constant (in fact lOOOM©) in 
our simulations (as Ale — 2500Mq in all our simulations). 

From its definition, /3 describes star formation efficiency 
in a certain sense. In fact, there is an explicit relation be- 
tween /3 and 77 (the commonly defined star formation effi- 
ciency), 



M,,,4o(r-|o + a')'/' \v - ^ 



(5) 



We find that our simulation results on the effect of the 
parent cloud on the cluster is better described by /3. Further 
discussion will be given later and the 40% Lagrangian radius 
of cluster will be denoted as r4o hereafter. 



2.1.3 Initial length scale ao 

The length scale of the Plummer potential describes the 
compactness of the cloud. The more compact the cloud is 
the more its effect on the cluster is when it disperses. We 



consider ao from 0.125 pc (compact) to 2.5 pc (loose) as 
listed in Table [1 

Each set (a,M6,ao) represents a model, and we name 
the model by the indices in Tables [T] [2] & |3] For example, 
BlOh corresponds to dispersion rate a — 0.2 (i.e., te — 1.5 
Myr), total cloud mass Alb ~ 5Mc = 125OOM0, and initial 
length scale of cloud ao = 1 pc. 



2.2 Model for the cluster 

2.2.1 Initial conditions 

The star cluster is prepared according to a Plummer distri- 
bution both in physical positions and velocities, which are 
required to achieve virial equilibrium. Note that the fact 
the cluster is in virial equilibrium does not mean that it i s 
also in dynamic al equilibrium (see e.g., iLada et all (|l984h . 
iGoodwinI (| 19971 )). In fact, right after we start the simula- 
tion, the cluster oscillates (shrinking and expanding) for a 
few times before settle down to a smoother (and more "natu- 
ral") distribution (see Fig. [3]). We, therefore, generate initial 
conditions according to the following steps: 



• generate a cluster with a Plummer distribution and a 
size about 1 pc; and we called this initial condition IC-0; 

• put the cluster into a molecular cloud, represented as 
a Plummer potential, with its centre of mass coincides with 
the centre of the potential; 

• turn off the cloud dispersion (i.e., a — 0), and run the 
code to two relaxation times; we called the conditions at one 
and two relaxation times IC-I and IC-II, respectively. 

To run the simulation, we turn on the cloud dispersion (i.e., 
a 7^ 0), after zero, one or two relaxation times accordingly. 

Fig. E] shows the ly-projection of the spatial and ve- 
locity distributions of different initial conditions. For IC-0, 
stars are restricted within 1 pc, and the velocity distribution 
has a very sharp upper limit. For IC-1 and IC-II, the cluster 
has smoother spatial and velocity distributions. 



2.2.2 Stellar mass function 

iGoodwinI (jl997h mentioned that there is no significant dif- 
ference between including a stellar mass function or not. 
However, for completeness and for comparison, we consider 
clusters with and without mass function. In both cases, the 
number of stars is 2500 and the total mass of the cluster 
is 2500 Mq. In cluster with equal mass stars, each star is 
1 Mq. In cluster with a stellar mass fu nction, w e adop t 
Salpeter mass function with slope -2.35 (jSalpeteil (|l955l ll 
and mass range from 0.32 to 32 Mq. 
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Table 4. Typical cases 
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Figure 3. xy-projcction of the spatial, mass and velocity distri- 
butions in radius of stars for the initial conditions IC-0, IC-I and 
IC-II, which correspond to preparing the cluster conditions at 0, 
1 and 2 relaxation times, respectively (see text for details). 
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Mb 
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3 RESULTS 

In order to investigate the problem as thoroughly as we can, 
for each dispersion rate a, each initial condition and each 
mass function, we perform four hundred simulation runs on 
the parameter space (Mi,,ao) as listed in Tables [2] & [S] We 
worked out five dispersion rates (see Table [1} , three initial 
conditions (see ii2.2.ip . and two mass functions equal mass 
and Salpeter mass funtion, see i]2.2.2|) . 

To illustrate the main results, we present four typical 
cases as listed in Tabled 

• B14s, Afb = 7.0Mc, ao = 2.375 pc (loose cloud); 

• B07h, Mh = 3.5Mc, ao = 1.0 pc (intermediate cloud); 

• BlOg, M(, = 5.0Mc, ao = 0.875 pc (intermediate cloud); 

• B18e, M(, = 9.0Afc, ao = 0.625 pc (massive and com- 
pact cloud). 

All of these cases have a = 0.2 (i.e., — 1.5 Myr), with 
initial condition IC-0 and Salpeter mass function. 

Fig.|4]shows how the half mass radii, Vh, vary with time 
for these four cases. In B14s (dash-dotted line), increases 
only a little bit. In B07h (dashed line), increases more 
than 2 times and becomes stable after 10 Myr. BlOg (dotted 
line) behaves similar to B07h, but the time to become stable 
is longer. Generally, increases in the first 5 to 10 Myr and 
then shrinks back. Cluster is considered survived in these 
three cases (B14s has a compact core, and B07h and BlOg 
are called loose). In BlSe (solid line), increases almost 
linearly. Cases such as BlSe, Vh will never decrease again and 
cluster is considered destroyed. More quantitative criteria 
for compact core-, loose- and destroyed-clusters will be given 
below. 



Figure 4. The evolution of half-mass radius r^. Solid line is the 
destroyed case BlSe. Dash-dotted line is the case with a compact 
core B14s. Dashed line and dotted line are loose cases B07h and 
BlOg, respectively. 

3.1 Effect of cloud dispersion 

The star cluster is in virial equilibrium in the potential well 
of the molecular cloud in the beginning. Every star is moving 
in a way according to how significant it is bound to the 
cloud and other stars. As the cloud is dispersing, its potential 
well is getting flatter and fiatter, and its effect is negligible 
after 10 Myr. During this process, the cluster expands. Some 
stars may escape. If the expansion is not too serious, the 
cluster could adjust and return to equilibrium after some 
time, otherwise, it is on the way to destruction. From the 
velocity distribution, we can describe the evolution of the 
cluster in four stages as shown in Fig. [5] In the figure, the 
dots are stars, the solid lines are tracks of elliptical orbits 
which semi-major axes are 10, 15, 20, 25 and 30 pc. 

(1) bound: stars orbiting the centre of mass of the cluster 
at high speed because of the additional gravity provided by 
the cloud; 

(2) expansion: during (and after) cloud dispersion, stars 
seem to escape; 

(3) inner part formation: stars with shorter periods may 
pass their aphelion and return (following the solid lines 
which describe the Keplerian orbits); 

(4) final structure at 30 Mry: stars with longer periods may 
return, and return stars will forget their previous orbits due 
to the complex interaction in the inner part. 
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Figure 5. DifTcrcnt evolutionary stages of case ElOi in a=0.5 
(te=0.625 Myr). Solid lines are tracks of elliptical orbit which 
semi-major axes are 10, 15, 20, 25 and 30 pc. There are four 
stages: (1) bound, (2) expansion, (3) inner part formation, and (4) 
final structure at 30 Myr. Red solid and dash lines in (4) seperate 
inner part, return stars and stars do not reach their aphelion. 



However, not every star could return if the dispersion of 
cloud is too serious. Our results show that the final structure 
can be grouped into three groups: (a) destroyed, (b) loose, 
and (c) compact core. Fig. [6] shows the three groups in the 
parameter space {Mf,,ao). By examining the spatial distri- 
bution, we find that the boundaries can be well described 
by a parameter e which we called the expansion ratio. It is 
defined as the ratio between the final and initial r4o of the 
cluster. 



r4o,i 



(6) 



The boundaries are e ~ 2 and 10. These numbers are chosen 
naively by examing the spatial distribution by eyes. 

(a) destroyed (e > 10): when clusters could not adjust 
themselves back to equilibrium in time, they will expand 
forever; e.g., case B18e; 

(b) loose (10 > e > 2): most clusters could adjust back to 
equilibrium but may become looser and loose part of stars; 
e.g., cases B07h and BlOg; 

(c) compact core (2 > e): in some cases, r4o expands 
slightly (less than twice), they become stable and develop 
a denser, e.g. case B14s. 

Fig. [7] shows the a;j/-projection of the final spatial distribu- 
tion of the cluster of the cases B18E, B07h, BlOg and B14s. 
Note that they have the same cloud dispersion rate. 
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Figure 6. Groups in the cases with te= 1.5 Myr and Salpeter 
mass function, (a) Three generic groups: destroyed (D), loose (L), 
and compact core (C) are shown in the parameter space (Afj,, ag). 
The dashed lines are the expansion ratio e = 2 and 10. The solid 
lines are cluster-cloud mass ratio f3 = 0.40 and 0.19, respectively. 
The two sets resemble each other very well (they almost overlap 
with each other), (b) Relation between e and /3. Each points on 
the figure is a simulation run (there are 400 points). It appears 
linear for e < 5, there is a lack of data in the interval 5 < e < 10 
and for e > 10 some other relation takes over. 
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Figure 7. xy-projection of the spatial distribution of the cluster 
at 30 Myr. There are three generic types: destroyed (B18e), loose 
(B07h & BlOg), compact core (B14s). 



3.2 Cluster-cloud mass ratio and expansion ratio 

at r-4(j 

As the cluster is bound by the parent cloud and its mem- 
ber stars initially. One would expect the more compact the 
cloud is the more the cluster will expand when the cloud 
is dispersed. After some experimentation, we find that the 
best parameters to delineate this anti-correlation is the ini- 
tial cloud-cluster mass ratio /3 and the expansion ratio e, 
both define at the 40% Lagrangian radius of the cluster r^o. 
Fig. [SJa shows how tight the correlation is. We should point 
out that if we define the cloud-cluster mass ratio and the 
expansion ratio with respect to other Lagrangian radii (say, 
30% to 50%), we get qualitatively similar results but the 
correlation is not as good as the 40% Lagrangian radius. We 
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Figure 8. Histogram of e for ie=l.l Myr. (a) Histogram from 
e=0 to 100. About 300 cases (out of a total 400) are under 10. 
(b) Histogram from e=0 to 10. Most are under 2 and the number 
drops off as e increases. 



Table 5. The value of f} corresponding to e = 2 and 10. 

te [Myr] 3.3 1.5 1.1 0.75 0.625 

02 (/3 at e = 2) 0.40 0.40 0.40 0.45 0.50 
/3lO (/3 at e = 2) 0.13 0.19 0.22 0.22 0.30 



still do not understand what is so special about r4o in our 
system. 

It is worthwhile to mention the statistics of e in our 
simulation sample (with same dispersion rate and mass func- 
tion) . A typical result is shown in Fig. [51 which shows the 
histogram of the number of runs resulting in e. Clearly, most 
of the runs lie in e < 10 and for e > 10 the data seems to be 
evenly distributed (Fig.|8K). When we zoom into < e < 10, 
we find that there is a drop off in 2 < e < 5, and within 
5 < e < 10 there is almost no data (Fig. [Sja). It seems that 
there is a transition region between 5 to 10. 

3.3 Dispersion rates 

We work out five dispersion rates of the cloud with e-fold 
time U is from 0.625 to 3.3 Myr (see Table[lJ. Fig.[9]shows 
the r4o expanding ratio e of three dispersion rates (with e- 
fold time = 3.3, 1.1, 0.625 Myr) with the same initial 
condition (IC-0). The e ~ 2 lines of the cases te — 3.3 and 
1.1 Myr almost overlap. As expected the region for destroyed 
cluster (e > 10) is larger for te=0.625 Myr (the most rapid 
dispersion rate in our simulations). In any case, many (more 
than two-third) of the clusters survive and remain intact 30 
Myr after the cloud started to disperse. 

Fig. [To] shows the relation between e and (3 for all five 
dispersion rates. Except te=0.625 (a=0.5), all of them look 
alike when e < 5 (dashed line is e = 5). Even for te— 0.625, 
the trend {/3 is larger when e is smaller) is the same as the 
others. 

From Fig. [To] we can read the corresponding value of /3 
for e = 2 and 10 for the five dispersion rates. The result is 
fisted in Table [5] 



3.4 Initial conditions in virial equilibrium 

Similar works usually start the cluster from virial equilib- 
rium state. However cluster in virial equilibrium might not 
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Figure 9. The three generic groups in different dispersion rate 
(cf, Fig.jeji). Solid line is for a=0.5 ( ie=0.625 Myr), dotted line 
is for a=0.3 ( te=l.l Myr) and dashed line is for o=0.1 (te=3.3 
Myr). For a=0.1 and 0.3, the e ~ 2 lines almost overlap, and 
although there is a separation in the e = 10 lines, they just look 
like a shifting. For a=0.5, the most rapid dispersion, although it 
looks a little different than the other two, wo can still divide it to 
the three generic groups. 
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Figure 10. Relation between e and /3 for different dispersion 
rates. The relations are fairly similar for the cases a = 0.1, 0.2, 
0.3 and 0.4. It is kind of different for a = 0.5, but we can still use 
the solid lines e = 2 and 10 to separate the destroyed, loose and 
compact core groups. 



necessary in dynamical equilibrium (e.g, iGoodwinI (|l997l ) 
and our tests, a brief description of our tests is given in 
i]2.2.1|l . We would like to know whether this would make a 
difference. Thus we consider three different initial conditions 
according to the prescription mentioned in i]2.2.1l (IC-0. IC- 
1, IC-2 for conditions taken at 0, 1, 2 relaxation times). 
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Figure 11. The relation between a and (3 (or te and /3) for e = 2 
(diamond and the upper line) and 10 (triangle and the lower line). 
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Figure 12. The three generic groups for te = 0.625 Myr in 
different initial conditions. Solid line, dotted line and dashed line 
are for IC-0, IC-I and IC-II, respectively. There is no significant 
difference between these three models. Nonetheless lC-0 has the 
smoothest boundaries. 



Fig. [12] shows the results of te — 0.625 Myr for diflFer- 
ent initial conditions. The results agree with each other rea- 
sonably well. Notheless IC-0 gives the smoothest divisions. 
Similar results are seen in other models. 



3.5 Stellar mass function and mass segregation 



Previous 
iBastian"^ 



works, such a s iGoodwin fc BastianI (|2006t ) 



GoodwinI (|2006l ). claimed that simulations with 



equal mass and with mass function should give similar out- 
come. For completeness, we re-examine this issue. All the 
results presented so far are cases with Salpeter mass func- 



Figure 13. The three generic groups in different mass functions. 
The dispersion c-fold time of the cloud is te = 1.5 Myr. Solid line 
is clusters with Salpeter mass function and dotted line is clusters 
with equal mass. The lines of e = 10 agree well with each other, 
but the lines of e = 2 are quite different, which is supposed to be 
a result of mass segregation. 



tion (see i )2.2.2p . We repeat all the simulations for clusters 
with equal mass (i.e., each star is 1 Mq). Note that the num- 
ber of stars and the total mass of the cluster are the same 
as those clusters with mass function (see ij2.2.2|l . 

The three generic groups: destroyed, loose and compact 
core are still applicable in the case of equal mass clusters. 
Fig.[T3]compares the results with and without mass function 
in the same dispersion rate and the same initial condition 
(IC-0). In the figure, the solid lines are clusters with Salpeter 
mass function and the dotted lines are clusters with equal 
mass. While the boundaries seperating the loose group and 
destroyed group for the two mass functions agree well with 
each other, the lines seperating the loose group and compact 
core group do not. This is true for every dispersion rate. 
We conclude that only the boundary between loose group 
and destroyed group does not depend on mass function. The 
boundary between loose group and compact core group are 
different and may be attributed to mass segregation. 

The concentration ratio for surviving clusters should be 
affected by mass function. Even for lower cluster-cloud mass 
ratio P could form a dense core in model with mass function. 
Examining mass distribution at later time, one should find 
that there is (and should be) mass segregation for cases with 
mEiss function. Fig. 1 141 shows the mass function of case C14s 
at 30 Myr (the end of our simulation). Since we do not have 
stellar evolution, there is no mass loss during the simulation. 
The mass function for the whole cluster remains Salpeter 
(solid line). In the inner region (1 pc, dashed line) the mass 
function is shallower, while in the outer region (between 1 to 
1.5 pc, dotted line) the mass function is steeper or at least 
the region has less massive stars. 

To learn how dispersion rates and initial conditions af- 
fect mass segregation, we plot the slope of the mass function 
against the cluster-cloud mass ratio /3 (Figs. [T5lfcll6|l. The 
general trend is the mass function is shallower when (3 is 
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Figure 14. Mass segregation of case C14s at 30 Myr (tiie end of 
our simulation). Solid line, dashed line, dotted line are the mass 
functions for the whole cluster (Salpeter), for stars within 1 pc 
(shallower than Salpeter), for stars between 1 pc to 1.5 pc (steeper 
than Salpeter), respectively. Note that all the first bins are scaled 
to the same number. 



Figure 15. Slope of mass function (within initial r^) against 
cluster-cloud mass ratio for different dispersion rates. The mass 
function of initial condition is Salpeter with slope -2.35. The gen- 
eral trend is more serious mass segregation when /3 is smaller. 
The scattering in small fj might be due to the number of stars is 
small at 30 Myr (the end of simulation). 



smaller, i.e., mass segregation is more serious at smaller (3. 
This can be understood as less massive stars are more likely 
to escape. (Recall that the cluster-cloud mass ratio /3 and 
the r4o expansion ratio e are anti-correlated, see Fig. 1101 or 

Eh.) 

In Fig. 1151 there is little difference between different dis- 
persion rates when j3 > 0.5. For smaller f} (i.e., when cloud 
dominates over cluster in the beginning), the data scatter a 
lot. The reason might be the number of stars remains in the 
cluster is small in these cases, and the statistics is not very 
good. 

In Fig. 1161 besides more serious mass segregation at 
small P, there is another increase in mass segregation when 
P ^ 0.8 for IC-I and IC-II. The reason can be traced back to 
the preparation of the initial conditions. For IC-I and IC-II, 
the cluster has already had some mass segregation in the 
beginning. It segregates more when /3 is large. After cloud 
dispersion, it becomes even more segregated. 



3.6 Infant mortality 

It is conceivable that the dispersion of a molecular cloud 
may destroy its embedded infant star clusters. The ques- 
tion is how effective is the process. From our simulation re- 
sults, we deem that the process is not very effective. Infant 
clusters can be destroyed only when the r^o cluster-cloud 
mass ratio P is below some value. Specifically, /3 < /3io (/3 
at e = 10) in Table [5] (also see Fig. Illf) . From the cases 
we considered, many infant clusters survive (some of them 
are loosened while some develop a compact core). Table [6] 
presents the frequency of the three groups (destroyed, loose, 
compact core) for different cloud dispersion rates and mass 
functions. Roughly 20% to 40% (e-fold time from 3.3 Myr 
to 0.625 Myr) will be destroyed. More than 30% for equal 
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Figure 16. Slope of mass function (within initial r^) against 
cluster-cloud mass ratio for the three initial conditions ( i|2.2.1l l. 
The cloud dispersion e-fold time is = 1-5. The general trend is 
more serious mass segregation when /3 is smaller. However, when 
(3 ^ 0.8 there is another increase in mass segregation for IC-I and 
IC-II (see text for a possible explanation). 

mass clusters and more than 40% for Salpeter mass function 
develop a compact core. 

We should point out that the boundaries {e — 2 and 10) 
between the three groups compact core, loose and destroyed 
are determined partly by eyes (from the spatial, velocity and 
mass distribution such as in Fig. [7]) and partly by numbers 
(such as the histograms in Fig. [Sjl . The boundaries between 
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Table 6. Infant mortality in % 



Salpeter 



te ]Myr] 


3.3 


1.5 


1.1 


0.75 


0.625 


Compact core 


55 


54 


54 


49.25 


42.75 


Loose 


26 


22 


16.75 


21.75 


21 


Destroyed 


19 


24 


29.25 


29 


36.25 


equal mass 












te [Myr] 


3.3 


1.5 


1.1 


0.75 


0.625 


Compact core 


37.5 


36.5 


35.5 


32.5 


30 


Loose 


45 


42.25 


38.5 


38 


37.75 


Destroyed 


17.5 


21.25 


26 


29.5 


32.25 



loose and destroyed (e = 10) we pick may not be precise. 
However, we claim that the exact location of the bound- 
ary (if there is an extra one) does not affect the estimated 
destroyed frequency listed in Table [S] As we discussed in 
i]3.2l there is a transition zone between e = 5 to 10 where 
only very few cases exist (see Fig. (Sj. As the boundary is 
somewhere around e « 10 (within the transition zone), the 
exact location will not matter too much on the estimated 
frequency for the two classes (destroyed and loose). 

We should point out that the range of the transition 
region is actually slightly different for different dispersion 
rates (see Fig. llOp . In most cases, the transition region runs 
from e = 5 to 15 (for slow dispersion, such as te = 3.3 Myr) 
and 5 to 30 (for rapid dispersion, such as te = 0.625 Myr). 
The exact location does not matter too much on the results 
as there are only very few cases inside the transition region. 



4 SUMMARY AND REMARKS 

Can a bound star cluster in a molecular cloud survive when 
the parent cloud is dispersed? We studied this problem by 
means of N-body simulations. The dispersion of the cloud 
makes the cluster to expand or even dissociate as the binding 
energy decreeises. However, our simulations show that this 
process is ineffective in destroying the bound cluster. 

The dispersion of the cloud is modelled by a Plummer 
potential with expanding length scale. We numerically simu- 
lated the behaviour of a bound cluster in this potential up to 
30 Myr. We performed numerous simulations and concluded 
that there are three groups of final morphology: compact 
core, loose and destroyed. It turns out that the 40% La- 
grangian radius expansion ratio of the cluster e (the ratio 
of r4o at 30 Myr to initial, see Eq. ([6])) is a good indicator 
for the final fate of the cluster. More importantly, we found 
that the initial r4o cluster-cloud mass ratio P (see Eq. (|4])) 
is tightly correlated with e. Hence we can "predict" the fate 
of the embedded bound cluster by examining its /3. 

The survival probabilities are high when f3 is high. Af- 
ter 30 Myr, more than 70% of the clusters still remain a 
shape of clusters, with a number density higher than 15/pc'^ 
within 1 pc radius. Different cloud dispersion rates provide 
similar results and even the largest rate in our simulations 
(e-fold time te = 0.625 Myr) does not disrupt all the clus- 
ters. Systems with and without mass function have different 
final densities but agree with each other well. 



We conclude that the infant mortality rate should be 
low if the embedded clusters are bound from the beginning, 
unless /3 is less than about 0.13 (for slow dispersion, such 
as te — 3.3 Myr) to 0.30 (for rapid dispersion, such as te = 
0.625 Myr) (see Table [Sjl. 

Near infrared observations indicate that the survival 
probability fo r emb edded clusters is about 4% to 7% 
iLada fc Ladal (|2003l ')). How should we reconcile our result 
with the observations? We deem that the dispersion of par- 
ent cloud can not account for the destruction of embedded 
bound clusters. Therefore, other more effective mechanisms 
must be responsible or the stars formed in m olecular cloud 
are not bound (see e.g.. iBonnell et al] (|2006l )). 
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